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FOREWORD 


This  report  represents  the  Second  Quarterly  Progress  Report,  1  September 
1962,  to  30  November  1962,  on  Contract  No.  AF  04(6ll)-8503,  a  contract  between 
Edwards  Air  Force  Base  and  The  University  of  Michigan.  The  aim  of  this  con¬ 
tract  is  to  investigate  the  feasibility  of  a  rotating  detonation  wave  rocket 
motor. 


Personnel  associated  with  the  various  phases  of  the  program  as  they  are 
divided  in  the  report  are  as  follows: 

I-A  -  J.  Fu 

B  -  M.  Sichel,  T.  David 
C  -  T.  C.  Adamson,  G.  Olsson 

II-A  -  G.  L.  Cosens,  J.  Brown 
B  -  G.  L.  Cosens,  K.  Ragland 
C  -  G.  L.  Cosens,  F.  Cheslak,  J.  Brown 
D  -  E.  Kurath,  S.  Schmidt 

This  project  is  directed  by  Professors  J.  A.  Nicholls  and  R.  E.  Cullen 
of  The  University  of  Michigan.  The  Air  Force  Project  Engineer  is  Richard 
Weiss  (DORR),  6593d  Test  Group  (Development),  Edwards  Air  Force  Base,  Cali¬ 
fornia. 
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SUMMARY 


This  report  presents  the  work  accomplished  on  the  rotating  detonation 
wave  engine  feasibility  program  during  the  period  1  September  1962  to  30  Novem¬ 
ber  1962. 

The  work  described  in  this  report  represents  the  continuation  of  the  ex¬ 
perimental  and  theoretical  studies  that  will  consume  approximately  the  first 
half  of  the  program.  The  last  half  of  the  program  will  be  devoted  to  the 
fabrication  and  testing  of  a  nominal  1000  lb  thrust  Ha-02  rocket  motor  utilizing 
the  detonative  mode  of  combustion. 

The  theoretical  studies  have  been  continued  in  the  two  areas  studied  during 
the  first  quarter.  In  addition,  the  study  of  a  third  problem  has  bean  initiated. 

(1)  The  study  of  the  structure  of  a  detonation  wave  through  a  two -phase 
medium  has  continued  emphasizing  the  droplet  shattering  mechanism.  Utilizing 
a  theoretical  relation  for  a  shear-type  breakup  developed  by  Dodd,  it  is  con¬ 
cluded  that  the  distance  predicted  for  droplet  breakup  behind  the  strong  nor¬ 
mal  shock  wave  associated  with  an  H2-O2  detonation  can  be  much  smaller  than  ex¬ 
perimentally  observed  distances  between  the  shock  and  the  zone  of  significant 
chemical  reaction.  It  is  concluded  therefore  that  the  droplet  shattering  effect 
occurring  behind  the  shock  wave  can  be  of  extreme  importance  in  substaining  a 
Chapman-Jouguet  detonation  in  dilute  sprays . 

(2)  The  study  of  the  heat  transfer  to  the  wall  of  a  rotating  detonation 
wave  engine  has  been  continued.  For  the  case  of  dissociative  recombination  at 
the  wall,  the  conclusions  are  that  the  heat  flux  to  the  wall  is  increased  about 
35  percent  over  the  case  without  dissociative  recombination  studied  earlier. 

Curves  are  presented  showing  this  effect  as  a  function  of  the  mole- fraction  of 
Ha  present  in  the  unburned  gases.  The  maximum  heat  flux  occurs  for  a  fuel-rich 
mixture  of  Xh2  =  0.75  (oxidizer  to  fuel  weight  ratio  of  5-5 )•  It  is  concluded 
also  that  according  to  the  simplified  model  derived  earlier,  that  the  heat  flux 
for  the  rotating  detonation  wave  engine  with  dissociative  recombination  occurring' 
still  is  of  the  same  order  as  the  heat  flux  at  the  throat  of  a  small  conventional' 
Ha-02  rocket  motor. 

(3)  The  study  of  an  analytical  model  for  the  gas  dynamics  associated  with 
the  rotating  detonation  wave  engine  has  been  initiated.  Utilizing  the  basic 
assumptions  that  the  flow  is  steady  with  respect  to  the  detonation  waves,  and 
that  the  flow  is  quasi -one  dimensional  in  the  circumferential  direction,  the 
conservation  equations  lead  in  general  to  a  system  of  first  order,  ordinary 
differential  equations  relating  the  flow  properties  in  the  chamber  and  a  di¬ 
mensionless  circumferential  distance,  q.  The  case  of  instantaneous,  complete 
mixing  of  the  burned  and  unburnsd  propellant  is  treated  as  well  as  the  case  of 
no  mixing. 
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Tha  experimental  studies  on  all  the  phases  reported  during  the  first  quarter 
are  continuing. 

(1)  Instrumentation  has  been  essentially  completed  on  the  100  lb  thrust 
motor  and  initial  runs  have  been  made .  Due  to  difficulties  with  the  miniature 
Kistler  pressure  transducer  for  recording  chamber  pressure  in  the  engine  as  a 
function  of  time,  the  transducer  has  been  installed  in  the  straight  detonation 
tube  for  shakedown  and  calibration  utilizing  a  known  wave  form; 

(2)  The  studies  of  the  pressure  and  temperature  effects  on  detonation 
velocities  have  progressed  to  the  point  of  substantiating  the  data  of  Moyle  down 
to  l60°K  at  a  pressure  of  one  atmosphere.  In  subsequent  tests,  it  appears  that 
the  coiled  tube  will  have  to  be  removed  from  the  cooling  bath  and  brought  up  to 
room  temperature  to  eliminate  the  ice  formation  of  the  combustion  products  on 
tube  walls  between  each  run. 

(3)  Instrumentation  of  the  condensation  chamber  for  the  study  of  droplet 
formation  by  rapid  expansion  and  condensation  is  almost  complete.  Tha  detona¬ 
tion  tube  has  been  fabricated  for  the  study  of  detonating  through  droplet  fields 
produced  by  either  condensation  or  injection  techniques. 

(4)  The  geometrical  tests  performed  to  date  have  investigated  the  wave 
structure  of  detonations  in  a  curved  tube  without  side  relief.  It  is  con¬ 
cluded  that  a  multiply-reflected  wave  occurs  in  the  burned  gases  behind  He -O2 
detonation  waves  in  nearly  stoichiometric  mixtures.  Also,  the  initial  wave 
front  is  slightly  inclined  to  the  radial  direction. 
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I.  THEORETICAL  STUDIES 


A.  DETONATION  WAVE  IN  HETEROGENEOUS,  LIQUID-GAS  MEDIA 

It  was  concluded  in  earlier  studies^  that  the  droplet  shattering  phenomenon 
could  be  of  considerable  importance  in  sustaining  stable  detonation  waves  through 
heterogeneous,  liquid-gas  media,  where  either  the  fuel  or  the  oxidizer  are  ini¬ 
tially  in  the  form  of  a  dilute  spray  of  liquid  droplets. 

The  following  is  a  study  utilizing  the  best  theoretical  and  experimental 
information  available  to  predict  the  order  of  the  breakup  times  for  relatively 
large  (500  -  1000  p)  liquid-oxygen  droplets  subjected  to  the  strong  shock  wave 
associated  with  a  H2-O2  detonation. 


1.  Some  Experimental  Observations  Concerning  the  Droplet  Breakup  Process 

4-6 

Experimental  results  show  that  two  types  of  droplet  breakup,  exist  in 
the  zone  behind  a  shock  wave  front:  the  bag -type  breakup  and  the  shear -type 
breakup.  The  successive  stages  of  a  droplet  breakup  process  can  be  described 
phenomenalogically:  (l)  For  a  bag-type  breakup,  the  increased  pressure  on  the 
droplet  surface,  due  to  the  high  convective  velocities  following  the  wave  front, 
flatten  the  droplet  into  a  disk-shape,  having  its  face  perpendicular  to  the  di¬ 
rection  of  the  flow.  Then  the  center  portion  is  blown  out  into  a  thin  hollow 
bag  anchored  to  a  heavy  rim.*  Finally  the  thin  bag  bursts  into  small  droplets, 
while  the  rim  breaks  into  fragments .  Photographic  evidence  on  this  phenomenon 
appears  in  References  4  and  5 •  (2)  For  a  shear-type  breakup,  in  contrast,  the 

droplet  is  distorted  into  a  saucer-shaped  disk  with  its  convex  surface  facing 
to  the  gas  flow.  Then  a  thin  ring  layer  is  drawn  out  from  the  edges  of  the 
disk.  Finally,  the  thin  layer  is  stripped  off  the  droplet  and  breaks  up  com¬ 
pletely  into  micro -droplets .  Engel^  showed  that  the  mechanism  of  micro-droplet 
formation  is  due  to  mechanical  effects  only  and  is  not  due  to  vaporization  of 
liquid.  Photographic  details  of  this  process  may  be  found  in  References  5,  6, 
and  7- 


2.  Various  Theories  on  Breakup  of  Droplets 

However  complex  the  phenomenon,  it  is  highly  desirable  to  treat  it  from 
both  the  experimental  and  theoretical  aspects .  The  following  is  a  discussion 
of  pertinent  theoretical  analyses  concerning  the  shattering  processes  of  a 
droplet : 


♦Lane 


5 


states  that  up  to  70  percent  of  the  liquid  remains  in  the  rim. 


1 


a)  Hlnze 1 s  Theory. —One  of  the  earliept  mathematical  analyses  on  this 
subject  was  presented  by  Hinze .®  The  theory  is  based  on  linearized  hydro- 
dynamical  equations  for  slight  deformations  of  a  droplet  in  an  air  flow.  The 
main  emphasis  was  in  the  derivation  of  the  relation  between  critical  speed  and 
critical  size.  Considering  the  influence  of  the  aerodynamic -pressure  on  a 
droplet  and  the  surface  tension  of  the  droplet,  Hinze  found  the  criterion  for 
determining  breakup  of  a  droplet  to  be  the  relevant  value  of  the  Weber  number. 

By  using  the  experimental  data  of  Merrigton  and  Richardson, ^  Hinze  indicated 
the  critical  Vfeber  number,  (We)cr  =  10.*  However,  the  critical  Weber  number 

can  be  derived  theoretically  by  using  the  critical  radius  (r)cr  [see  Eq.  (3)  of 
Section  b]  and  the  definition  of  the  Weber  number.  It  shows  the  critical  Weber 
number,  (We )cr  =  8,  which  agrees  closely  with  Hinze's  value.  Hinze  also 

found  that  the  effect  of  viscosity  of  liquid  on  deformation  appears  small  ex¬ 
cept  for  large  values  of  viscosity,  e.g.,  viscosities  of  the  order  of  glycerol. 

In  this  case,  the  viscosity  effect  tends  to  retard  the  breakup.  The  effect  of 
viscosity  on  breakup  has  been  studied  in  more  detail  by  Hanson.11  He  found 
that  when  the  viscosity  of  a  liquid  is  less  than  about  10  centistokes,  the 
effect  on  droplet  breakup  is  negligible;  while  if  it  is  between  10  and  100 
centistokes,  the  breakup  process  is  retarded. 

b )  Mathematical  Model  for  Bag  Type  Breakup  of  a  Droplet.  This  mathematical 
analysis  was  proposed  by  Gorden.12  in  the  analysis,  it  1b  assumed  that  a  cylin¬ 
drical  plug  of  diameter  r  and  length  2r  (where  r  is  the  radius  of  the  droplet) 

is  extruded  from  a  droplet;  while  the  rest  of  the  droplet  around  the  cylinder 
remains  at  rest.  This  deformation  is  caused  by  the  air  pressure  in  front  of  the 
droplet,  but  is  retarded  by  surface  tension,  viscosity,  and  inertial  forces  of 
the  droplet.  Estimating  the  magnitude  of  the  pressures  caused  by  these  forces, 
and  combining  them  with  the  inertial  effect,  a  differential  equation  for  the 
extruded  cylinder  is  obtained,  i.e.. 


dv  _1_/  1  ~  16  v  8S\ 

dt  =  p jD  \2  P8  "  D  -  Dy 


(1) 


where  v  =  velocity  of  the  ext. hided  cylinder.  Equation  (l)  is  solved  for  the 
instantaneous  velocity,  v,  and  the  resulting  equation  is  in  turn  solved  for  the 
instantaneous  displacement  of  the  cylinder  as  a  function  of  time.  The  displace¬ 
ment  is  then  set  equal  to  the  droplet  diameter  D  to  determine  the  total  breakup 
time,  tf  The  result  is: 


2(1 6uj)2 

P/D"(p/  '  ^  ) 


i6n/th 

• 1+exp 


(2) 


*Hinze®  defined  We  =  (pgl^rJ/S  but  Williams10  defined  We  *  [pgUe(2r )  ]/S. 
In  this  report,  Hinze's  definition  will  be  adapted. 
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where  t^  is  the  breakup  time .  The  breakup  time  approaches  infinity  if  a  droplet 
is  stable.  Hence  it  is  obvious  from  Equation  (2)  that  the  critical  condition 
occurs  at 


(V2 


l6s\ 

D  ) 


0 


which  implies  the  critical  radius  is  given  by 


(r) 


cr 


8s 

p  IP 

Mg 


(3) 


Equation  (5)  is  formally  identical  with  ^ane 's  empirical  formula,  which  has  been 
used  in  correlating  experimental  results^'^  between  critical  speed  and  critical- 
size.  However,  Equation  (3)  will  be  used  in  this  report  for  the  same  purpose. 


Since  Equation  (2)  cannot  be  solved  analytically  for  the  breakup  time,  a 
useful  approximate  solution  is 


(2-a) 


Hanson^1  and  co-worker  showed  experimentally  that  the  effect  on  breakup  can 
he  neglected  when  the  viscosity,  p/  <L  10  centistokes.  For  this  case  the  bag-type 
breakup  time  is  approximately 


% 


2D  p^l 2 


16s  N 1/2 


(2-b) 


It  should  be  noted  that  according  to  this  theory,  Eq.  (2-b)  can  predict  the 
breakup  time  of  droplets  larger  than  critical  size,  and  for  droplets  of  critical 
3izes  the  breakup  time  is  infinite. 

q)  Mathematical  Model  for  Shear -Type  Breakup  of  a  Droplet. — This  mathe¬ 
matical  model  is  attributed  to  Dodd. 13  The  basic  considerations  of  this  theory 
are  as  follows:  The  tangential  frictional  aerodynamic  forces  cause  an  internal 
circulation  in  the  surface  layer  of  a  droplet.  In  turn,  the  centrifugal  effect 
resulting  from  the  circulation  causes  bulges  to  form  on  the  surface .  On  the 
other  hand,  the  surface  tension  of  the  droplet  has  a  retarding  effect  on  bulge 
formation.  If  the  pressure  due  to  the  centrifugal  effect  equals  a  fraction,  F,  of 
the  surface  tension  pressure,  then  breakup  occurs.  The  assumptions  used  in  Dodd's 


T 
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analysis  are:  1.  The  circulation  velocity  is  uniform  across  the  thin  surface 
layer;  2.  Inside  the  moving  layer  the  liquid  is  at  rest;  3.  The  frictional 
force  on  the  surface  of  the  droplet  is  equal  to  O.332  ■\j(U3UgPg)/r  and  this 
force  acts  over  a  region  of  length  r  (see  Fig.  2  of  Ref.  l);  4.  The  thickness 
of  the  moving  layer  is  "fr";  where  0  <  f  <  1.  Based  upon  the  balance  of  the 
energy  transfer  rate  between  the  gas  flow  behind  the  shock  front  and  the  droplet, 
a  breakup  time  for  a  droplet  is  derived: 


tb  = 


2r 


(1.5U )3/a  °-532 


sp/  JF 

pgpg 


or 
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Sp! 

^gpg 


By  using  Engel's^  experimental  data,  with  \jF/f  *  0.5 • 


*b 


1.64r  SPj 

U3/2  <J  PgPg 


(4) 


Due  to  the  inequality  of  pressure  on  the  upstream  and  downstream  sides,  each 
bulge  is  subjected  to  an  aerodynamic  drag  tending  to  move  the  bulge  off  the 
surface  of  the  droplet.  Dodd1*  has  pointed  out  that  the  effect  of  the  aero¬ 
dynamic  force  which  reinforces  the  frictional  force  on  the  droplet  surface, 
will  decrease  the  breakup  time.  Therefore,  the  breakup  time  given  by  Eq..  (4)  ? 
is  too  high,  and  some  modification  factor  is  needed  to  give  the  correct  breakup 
time .  Inasmuch  as  Dodd  did  not  attempt  to  evaluate  this  factor,  an  attempt 
will  be  made  to  estimate  this  factor  approximately  by  comparing  the  order  of 
magnitude  of  the  aerodynamic  force  and  the  frictional  force . 

Let  p  =  the  aerodynamic  pressure  =  l/2  pgU2 

pf  =  the  frictional  force  per  unit  area  =  0-352  , (U3  PgPg)/r 


+  p  =  1  +  JL. 

Pf  Pf 
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Since  the  impulse  on  the  droplet  due  to  the  gas  is  equal  to  product  of  the  force 
and  the  time  over  which  it  acts,  it  implies  the  breakup  time  is  inversely  pro¬ 
portional  to  the  force.  Hence  breakup  time  given  by  Eq.  (4)  should  be  modified 
by  a  factor  l/(l  +  1-5  .JRe),  i.e., 


tb 


1.64r  SP  / _ 1 

U3/2  VgPg  l  +  1.5  Jte 


(4-a) 


In  investigating  the  disintegration  of  water  drops  in  an  air  stream,  where 
the  relative  velocity  between  air  and  the  drop  increases  gradually,  Dodd  found 
the  critical  radius  for  a  droplet  is 


(r)cr 


6.67S 

pgPua 


where  p  is  a  correction  factor  0  <  p  <  1.  Lane,  using  Taylor’s  loading  spring 
concept  suggests  that: 

Velocity  for  bursting  by  suddenly  applied  blast  _  0  71 

Velocity  of  bursting  in  steady  stream 

In  Reference  14,  Lane  found  this  factor  agrees  reasonably  well  with  experimental 
data.  Using  this  value,  the  correlation  between  critical  relative  velocity  and 
critical  size  becomes 


(r)cr  =  13S/ ( PgU2 ) 


(5) 


3.  Comparison  Between  Theories  and  Experiments 

In  order  to  estimate  the  validity  of  these  mathematical  models,  some  calcu¬ 
lations  have  been  made.  The  results  are  listed  in  Tables  I,  II,  III,  and  IV. 

Since  it  has  been  pointed  out  by  Rabin  that  small  droplets  undergo  a  bag-type 
breakup  while  the  large  droplets  experience  a  shear-type  breakup,  Eq.  (3)  is 
used  to  correlate  the  critical  conditions  for  small  droplets,  while  Eq.  (5)  is 
used  for  large  droplets.  When  applying  Gorden's  correlation,  Eq.  (3),  to  the 
experimental  data,  in  Reference  .6,  for  (r)cr  =  500  p,  the  deviation  between 
theoretical  and  experimental  critical  radii  is  over  50  percent.  It  is  felt 
that  Eq.  (3)  is  valid  only  for  r  <  500p,  while  Eq.  (5)  applied  for  r  >  500  u- 

Equation  (4-a)  gives  the  breakup  time  for  the  shear-type  breakup  for  droplets 
of  critical  size  as  well  as  for  droplets  larger  than  critical  size .  Some  calcu¬ 
lated  results,  by  using  Equation  (4-a)  and  the  experimental  data  given  by  Reference 
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(6)  are  listed  in  Table  III .  It  indicates  the  breakup  time  for  droplets  of 
r  r  500  4  and  relative  air  flow  velocities  of  the  range  55  ft/sec  <  U  <  93  ft/sec, 
is  about  1  millisecond.  The  velocities  used  in  the  above  experiments  were  much 
lower  in  comparison  to  the  velocity  occurring  behind  the  strong  normal  shock 
wave  associated  with  the  detonation  process.  Hence,  the  breakup  time  in  the  zone 
behind  the  detonation  wave  front  would  be  expected  to  be  much  smaller  than  the 
values  listed  in  Table  III.  In  order  to  estimate  the  order  of  the  droplet  break¬ 
up  time  in  the  detonation  case,  the  following  analysis  is  made  utilizing  Equation 
(4-a )  and  the  detonation  characteristics  of  hydrogen-oxygen  mixtures  calculated 
by  Moyle. 15  For  a  mixture  having  a  mole  fraction  of  O.667  hydrogen,  a  detonation 
velocity  of  about  9^00  ft/sec  is  obtained.  Assuming  that  the  structure  of  a  spray 
detonation  consists  of  a  normal  shock  wave  followed  by  combustion,  the  gas  veloc¬ 
ity  immediately  behind  the  normal  shock  front  is  supersonic  and  causes  a  detached 
shock  wave  in  front  of  the  droplet.  It  is  the  gas  flow  behind  the  detached  shock 
that  is  responsible  for  droplet  breakup.  Using  the  conservative  minimum  value 
of  the  sub -sonic  velocity  behind  the  detached  shock,  the  breakup  time  for  oxygen 
droplets  having  a  size  range  of  50°  ^  r  <L  1000  u  have  been  calculated.  The  re¬ 
sults  are  shown  in  Table  IV.  The  breakup  times  are  in  order  of  0.1  u-sec .  In 
this  time  interval  the  detonation  wave  front  will  travel  a  distance  approximately 
equal  to  1  x  10 *2  in.  This  distance  will  be  denoted  by  Lb-  Morrison^®  has  found 
from  schlieren  photographs  utilizing  fuel-lean  gaseous  mixtures  that  there  can 
be  a  measurable  distance  between  the  leading  shock  wave  front  and  the  zone  where 
significant  chemical  reaction  takes  place.  Under  these  conditions,  detonation 
waves  propagating  at  velocities  in  good  agreement  with  velocities  calculated  for 
Chapman-Jouguet  detonation  waves  give  evidence  that  this  zone  may  be  as  great  as 
two  tube  diameters  or  more  in  a  0.5  in.  diameter  tube.  Denoting  this  distance 
by  Lc,  it  follows  that  for  H 2-02  detonations  in  lean  gaseous  mixtures  Lq  =  0(l) 
and  therefore  Lb  «  Lc • 


4 .  Conclusions 

The  Weber  number  appears  to  be  an  important  criterion  concerning  the 
droplet -shattering  phenomenon.  Hinze  has  shown  there  exists  a  critical  Weber 
number  of  a  dilute  spray,  (We)Cr  =  10-  When  the  Weber  number  is  larger  than 
this  value,  the  droplets  tend  to  break  up. 

Utilizing  a  theoretical  relation  developed  by  Dodd,  the  breakup  times 
for  oxygen  droplets  in  a  H2-O2  spray  detonation  have  been  calculated.  For 
oxygen  droplets  within  a  size  range  of  500  ^  r  £  1000  n,  the  breakup  time  is 
of  the  order  of  0.1  nsec .  In  this  time  interval  the  detonation  wave  front  will 
travel  a  length.  Lb  =  1  x  10-2  in.  In  comparing  this  order  of  LB,  with  the 
order  of  Lc,  the  length  between  the  shock  wave  front  and  the  zone  of  significant 
chemical  reaction,  it  follows  that  Lb  <  <  Lc.  Thus  it  is  indicated  that  large 
droplets  will  be  shattered  into  microdroplets  in  a  zone  sufficiently  small  be¬ 
hind  the  shock  wave  front,  that  if  the  subsequent  evaporation  and  combustion  of 
these  microdroplets  occur  rapidly  enough,  a  heat  release  sufficient  to  sustain 
a  Chapman-Jouguet  detonation  wave  is  possible .  Thus  it  appears  that  the  droplet 
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breakup  processes  occurring  behind  the  shock  wave  can  be  an  important  mechanism 
in  stabilizing  a  detonation  wave  in  a  dilute  spray. 

It  is  obvious  that  experimental  results  substantiating  these  theoretically 
predicted  droplet  breakup  times  behind  strong  shock  waves  is  highly  desirable. 


B.  STUDY  OF  HEAT  TRANSFER  IN  THE  ROTATING  DETONATION  WAVE  ENGINE 


The  work  described  below  is  a  continuation  of  the  heat  transfer  study, 
which  was  presented  in  detail  in  Reference  1.  The  effects  of  dissociation  and 
recombination  and  of  fuel-oxidizer  mixture  ratio  upon  the  heat  transfer  have 
been  investigated.  As  mentioned  in  Reference  1,  the  recombination  of  dissociated 
species  which  are  present  behind  the  detonation  can  result  in  an  appreciable  in¬ 
crease  in  the  heat  flux  across  the  boundary  layer. 

In  a  dissociating  gas,  erasrgy  is  transported  not  only  by  conduction  but 
also  by  the  diffusion  of  atomic  species .  The  heat  flux  across  a  boundary  layer 
in  a  dissociated  gas  may  be  increased  by  the  recombination  of  atomic  species 
within  the  boundary  layer  and  by  diffusion  of  atomic  species  from  the  free  stream 
to  the  cold  wall  where,  in  most  cases,  recombination  with  heat  release  occurs. 

When  the  Lewis  number,  Le  =  (CpDp)/k,  is  unity  the  heat  transfer  across  a  boundary 
layer  in  contact  with  a  cold,  catalytic  wall  is  proportional  to  the  difference 
between  the  free  stream  and  wall  enthalpies  regardless  of  the  energy  transfer 
mechanism. 17  The  enthalpy  of  course  must  include  the  chemical  heat  of  formation. 

For  a  preliminary  estimate  of  the  increase  in  heat  transfer  due  to  dissocia¬ 
tion  it  has  been  assumed  that  Le  =1.0.  Consequently  it  has  been  possible  to  use 
the  boundary  layer  heat  flux  results  of  Reference  1,  the  only  change  being  in 
the  calculation  of  the  enthalpy  difference  across  the  boundary  layer. 


The  heat  flux  to  the  wall  also  depends  upon  the  ratio  of  hydrogen  to  oxygen 
in  the  explosive  mixture  into  which  the  detonation  propagates .  Consequently, 
calculations  of  wall  heat  flux  have  been  made  for  0.40  <L  Xh2  0 .80  where  Xjjg 
is  the  mole  fraction  of  molecular  hydrogen.  The  composition  of  the  combustion 
products  at  the  Chapman-Jouguet  plane  have  been  computed  by  Moyle, ^-5  and  his 
results,  which  were  used  in  the  calculations,  are  reproduced  in  Fig.  1.' 

The  heat  flux  was  computed  as  described  in  Reference  1,  i.e. 


(hr  -  hy)Ty 

(VD  -  Ve)Prm2/3 


(1) 


where  tw,  the  shear  stress  at  the  wall  is  given  by  Equation  (9)  in  Reference  1* 
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and  where  the  subscript  m  indicates  that  p  and  p  are  calculated  at  the  tempera¬ 
ture  corresponding  to  the  Eckert  reference  enthalpy 


hm  =  0.5  (hv  +  hg)  +  0.22  (hr  -  1^) 


(2) 


The  free  stream  enthalpy  per  unit  mass  of  gas  consisting  of  the  constituents 
of  the  equilibrium  mixture  at  the  Chapman -Jouguet  plane  was  computed  according 
to  the  equation 


-  -«1,298.16-K*4V 


i  =  1 


where  W  i8  the  average  molecular  weight  of  the  mixture  given  by 

V  -  £  XjWi 


(3) 


(4) 


the  enthalpies  and  standard  heats  of  formation  have  been  taken  from  the  tables 
in  Benner. The  recovery  enthalpy  of  the  turbulent  boundary  layer  was  once 
again  computed  according  to 


In  computing  the  reference  temperature  from  hm  it  was  assumed  that  the  fluid 
composition  was  the  same  as  in  the  free  stream.  The  validity  of  this  assump¬ 
tion  is  somewhat  questionable  and  must  be  studied  further. 

A  wall  temperature  of  9006  R  was  assumed.  Since  dissociation  is  negligible 
at  this  temperature  the  gas  at  the  wall  will  consist  of  He  and  HeO  for  the  rich  mix¬ 
tures  and  0g and H2O  for  the  lean  mixtures. 

The  results  of  the  calculations  described  above  are  presented  in  Table  V 
and  in  Fig.  2.  It  can  be  seen  that  the  heat  flux  values  computed  by  taking 
dissociation  into  account  are  about  55^  more  than  the  non-dissociative  value. 

The  maximum  heat  flux  occurs  for  a  fuel-rich  mixture  with  XRa  *  0.75  or  a  ratio 
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of  oxidant  to  fuel  weight  of  5.3.  Again  assuming  as  in  Reference  1  that  the 
heat  pulse  lasts  for  only  one  third  of  the  period  between  the  passage  of 
successive  waves,  the  maximum  value  of  the  overall  heat  flux,  Qy,  to  the  walls 
of  the  100-lb  thrust  engine  at  Xjj2  =0.75  becomes 


max 


1660  Btu/ft2  sec 


11.5  Btu/in.2  sec 


For  the  mixture  with  Xg2  =  0  60,  which  was  considered  in  Reference  1  the  overall 
heat  flux  density  would  be 


Qw 


XH2 


0.60 


1480  Btu/ft2  sec  =  10.2  Btu/in.2  sec 


It  is  interesting  to  compare  the  above  values  of  Q  to  the  nozzle  throat  heat 
flux  of  9-95  Btu/in.2  sec  computed  by  Curran  et  al^9  for  a  conventional  H2-O2, 
1000-lb  thrust  rocket  for  a  mixture  with  Xjj2  =  0.84  (corresponding  to  maximum 
specific  impulse).  Thus  the  heat  flux  of  the  100-lb  rotating  detonation  wave 
engine  is  of  the  same  order  as  the  throat  heat  flux  of  a  small  conventional 
engine,  and  this  is  still  in  essential  agreement  with  the  conclusion  reached  in 
Reference  1. 

Nevertheless,  the  calculations  above  indicate  basic  differences  between  the 
rotating  detonation  wave  and  conventional  rocket  motors .  For  a  given  mixture 
ratio  the  temperature  behind  a  detonation  is  higher  than  in  a  conventional  com¬ 
bustion  chamber  with  a  corresponding  increase  in  the  dissociation  of  the  combus¬ 
tion  products.  This  can  be  seen  from  Fig.  3  which  shows  the  variation  of  Te, 
the  temperature  behind  the  detonation  and  Tc,  the  temperature  in  a  conventional 
300  psia  combustion  chamber^0  with  Xj|2 .  At  the  mixture  ratio  for  which  Isp,  the 
specific  impulse  of  a  conventional  engine  is  a  maximum,  Tc  is  about  1000°F  lower 
than  Te.  The  dissociation  behind  the  detonation  is  thus  sufficient  to  cause  an 
appreciable  increase  in  the  heat  flux,  as  is  shown  by  the  calculations  above. 
Even  though  the  temperature  Te,  reaches  a  maximum  value  near  the  stoichiometric 
composition  with  Xh2  =  0.667,  the  velocity  difference,  (V])-Ve),  across  the  wave 
continues  to  increase  as  Xg2  increases .  Since  heat  flux  increases  with  both 
(VD-Ve),  and  (hr  -  hw),  the  maximum  heat  flux  occurs  with  a  rich  mixture, 

XHs  =  0  •75,  rather  than  occurring  near  the  stoichiometric  mixture  ratio, 

Xgg  =  0.667. 
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C.  ESTABLISHMENT  OF  AN  ANALYTICAL  MOEEL  OF  THE  ROTATING  DETONATION  WAVE  ENGINE* 


1.  Introduction 

It  is  highly  desirable  to  he  able  to  establish  the  effects,  of  the  various 
design  parameters  (engine  sizing,  flow  rate  of  propellants,  etc.)  on  the  theo¬ 
retical  performance  of  the  rotating  detonation  wave  engine.  To  accomplish  this 
it  is  necessary  to  formulate  an  analytical  model  of  the  engine.  This  model, 
while  retaining  the  essential  physical  features,  should  be  simple  enough  to  per¬ 
mit  insight  into  the  gasdynamics  of  engine  operation.  Also,  the  results  of  the 
analysis  should  correlate  with  data  from  the  experimental  engine  tests.  Once  a 
reasonable  analytical  model  is  established  it  can  be  used  as  a  design  tool  and 
as  a  basis  for  comparing  the  theoretical  performance  of  the  rotating  detonation 
wave  engine  with  that  of  other  types  of  propulsion  devices. 

The  fluid  flow  field  in  the  engine  is  three  dimensional  and  unsteady  with 
heat  and  mass  addition  and  turbulent  mixing  taking  place.  As  a  result,  simpli¬ 
fying  assumptions  must  be  made  to  achieve  a  solvable  mathematical  model.  The 
key  assumption  is  that  the  detonation  wave  velocity  is  constant,  so  that  the 
fluid  flow  field  is  steady  and  quasi-one  dimensional  with  respect  to  the  detona¬ 
tion  waves.  Both  the  burned  and  unburned  propellants  sure  treated  as  ideal  gases 
and  the  effects  of  shear  forces  and  heat  transfer  on  the  flow  are  neglected. 
These  restrictions  provide  the  basis  for  a  simplified  analytical  model  of  the 
engine.  The  effects  due  to  the  turbulent  mixing  of  the  burned  and  unburned 
propellants  are  to  be  estimated  by  considering  alternatively  the  limiting  cases 
of  instantaneous  complete  mixing  and  of  no  mixing. 


2.  Fundamental  Assumptions  and  Conditions 

Because  of  the  complex  nature  of  the  gasdynamics  of  the  rotating  wave 
engine  it  is  necessary  to  impose  restrictive  assumptions  to  establish  a  rea¬ 
sonable  mathematical  model.  The  following  is  a  presentation  of  these  assump¬ 
tions  with  remarks  concerning  their  validity. 

(a)  The  fluid  flow  field  in  the  engine  is  assumed  to  be  steady  and  quasi- 
one  dimensional  with  respect  to  the  moving  detonation  waves. 

The  necessary  conditions  implied  in  this  assumption  are: 

(1)  The  detonation  wave  speed  has  reached  a  steady  state,  constant  value. 

(2)  The  chamber  dimensions  are  small  compared  to  the  radius  of  curvature 


•This  analysis  is  largely  the  effort  of  T.  C.  Adamson,  Professor  of  Aeronautical 
and  Astronautical  Engineering  at  The  University  of  Michigan.  The  authors  are 
extremely  grateful:  for  this  contribution. 
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of  the  chamber  centerline . 


(3)  The  chamber  fluid  properties  vary  only  in  the  circumferential  direction. 

(4)  In  the  chamber  the  velocity  in  the  axial  direction  is  small  compared 
to  the  velocity  in  the  radial  direction. 

Order  of  magnitude  condition  for  the  radial  pressure  variation  in  the  chamber 
to  be  small: 


For  a  curved  streamline 


dP 

dn 


In  order  to  have  a  negligible  pressure  variation  in  the  radial  direction. 


Hence, 


But, 


^  «Z. 

dn  lc 


P 

pva2 


P  =  1 

PV  ~  ^ 


which  is  of  order  one.  Then  the  condition  for  the  radial  pressure  variation 
to  be  small  is 


*c 

R 


«  1 


Order  of  magnitude  condition  for  the  unsteady  fluctuation  of  pressure  in 
the  chamber  to  be  small: 


The  time  for  an  acoustic  wave  to  cross  the  chamber  in  the  axial  direction 
must  be  much  less  than  the  time  the  detonation  wave  takes  to  travel  the  distance 
between  successive  waves.  Therefore,  the  condition  is, 


a 


2nR 

NVW 
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or 


x„ 

-2.  « 


2na 

NVv 


The  minimum  value  of  a  is  ai  and  since 


Vw 

ai 


Mw 


The  condition  becomes 

2“.  «  S*. 

R  NM* 

Mw  is  of  the  same  order  as  2it  and  N  is  of  the  order  one .  Then  the  condition 
for  negligible  axial  unsteady  fluctuation  of  pressure  in  the  axial  direction  is 


For  a  given  engine  geometry  these  comparisons  can  be  used  to  check  the 
assumption  of  one  dimensional  flow. 

(b)  The  fluid  flow  through  the  convergent  section  of  the  nozzle  is 
assumed  to  be  quasi -one  dimensional  in  the  axial  direction. 

As  a  consequence,  the  axial  component  of  velocity  at  the  throat  is  equated 
to  the  local  speed  of  sound  and  the  circumferential  velocity  component  is  assumed 
to  be  invariant  with  respect  to  the  axial  direction. 

(c)  It  is  assumed  that  the  detonation  wave  is  a  plane  discontinuity. 

The  detailed  structure  of  the  detonation  wave  is  not  to  be  considered.  It 
is  to  be  treated  hydrodynamically  as  a  shock  wave  with  heat  addition. 

(d)  Both  the  burned  and  unburned  propellants  are  assumed  to  be  ideal  gases. 

The  influences  of  the  variation  of  gas  parameters  with  temperature  and 
pressure  are  to  be  neglected  and  all  chemical  reaction  is  assumed  to  take  place 
only  at  the  detonation  wave . 

(e)  Shear  forces  and  heat  transfer  effects  are  to  be  neglected. 

(f )  The  propellant  mass  flow  rate  through  the  injector  is  assumed  to  be 
constant  and  uniform  in  distribution. 
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Assuming  the  propellants  are  delivered  to  the  injector  in  the  gaseous  state, 
this  condition  is  satisfied  if  the  injector  to  chamber  pressure  ratio  is  always 
at  or  above  the  critical  value  and  if  the  injector  design  is  axially  symmetric. 

In  order  to  derive  the  conservation  equations  (continuity,  momentum,  and 
energy)  two  special  cases  are  designated.  In  the  first  case  the  burned  and  un¬ 
burned  gases  are  treated  as  being  completely  mixed,  and  in  the  second  case  the 
gases  are  treated  as  being  unmixed.  It  is  believed  that  this  approach  will  give 
some  insight  into  the  effect  of  mixing  on  the  performance  of  the  engine. 


3.  Definition  of  Coordinate  Systems 


Because  the.  radial  width  is  considered  to  be  small  compared  to  the  radius 
of  curvature,  the  chamber  can  be  "unrolled"  and  stretched  out.  Then  the  problem 
is  one  wherein  a  long  tube  of  constant  cross-sectional  area  has  unburned  mass 
entering  on  one  side  and  burned  mass  being  ejected  on  the  other.  In  this  tube 
detonation  waves  occur  periodically  a  distance  L  apart,  each  moving  at  a  con¬ 
stant  velocity,  Vw.  From  the  viewpoint  of  an  observer  on  any  wave  the  problem 
becomes  steady  and  periodic. 


The  Wall  Fixed  System: 


x 


Note :  Arrows  define  the  direction  of  positive  velocity. 
The  two  systems  are  related  by  the  transformations 
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v  =  Vw  -  va 


y  =  vw t 

where  t  *  0  immediately  after  a  wave  passes. 

4.  Jump  Conditions  Across  the  Detonation  Waves 

The  fluid  flow  properties  at  y  =  0  and  y  =  L  are  related  by  the  hydrodynamic 
Jump  conditions  across  the  detonation  wave : 


PlVl  =  p0  vq 

Pi  +  Pivi2  =  PD  +  P0  vo2 
CPl  Tx  +  i  Vi2  +  nxQ  *  Tn  +  i  v„2 


Po  xo 


are . 


With  the  ideal  gas  assumption  the  thermal  and  caloric  equations  of  state 


_  Rq 

P  =  p  —  T 
m 


7  R0T 

C-nT  =  — Z - 2- 

p  7  -  1  m 


and  the  speed  of  sound  is  given  by 


a2  .  2S. 

P 


The  Jump  conditions  can  be  expressed  in  the  non-dimensional  form 

PlVl  =  1 

Pi  +  70pivi2m02  =  l  +  70Mq2 


fPl  ^ - 2  u "  -  7°'1  5 


Vi2  Mq  +  (i!§  =  1  + 


2  Mo 


mo  _  - 
Pi  =  ="  Pi  Ti 

mi 


CpiT,  = 


7o 


-  1 


ai2 


7i  -  1 

u A 

7o  Pi 


Ik 


These  equations  may  be  solved  to  give  Pi,vi,  etc.  in  terms  of  MQ  and  Q. 
For  example , 


^i  *  b  +  Jb2  +  c 


where  b  and  c  are  defined  by 


b  =  LL 
To 


LUs^sl. 


(n  +  i)Mo£ 


-g  (lx  -  ± 

=  r0  -  i  \7i  +  i>  m0- 


mQ  -  Z2~^1  Mq2  -  i| 


If  Mq  =  1  and  >1  =  yQ  =  y  this  becomes 


Vl 


1  + 


fgifiS 

sir  +  i 


In  terms  of  Vi, 

Pi  =  1  -  70M02(v!  -  1) 

ai2  =  7-±  [1  -  70Mo2(vi  -  l)]vi 

y  o 


5  .  Equations  for  the  Case  of  Complete  Mixing  of  the  Burned  and  Uriburned 
Propellant 

In  this  case  it  is  assumed  that  the  burned  and  unburned  gases  mix  completely 
and  instantaneously  and  that  the  specific  heat  ratios  and  the  molecular  weights 
have  one  average  value . 

Choose  a  small  increment  of  tube ,  of  length  dy,  and  write  the  conservation 
equations  for  two  control  volumes:  one  bounded  by  the  chamber  and  the  second 
bounded  by  the  convergent  section  of  the  nozzle.  In  deriving  the  functional 
relations  for  the  chamber  conditions  it  is  not  necessary  to  consider  the  diver¬ 
gent  section  of  the  nozzle  .  This  aspect  is  to'  be  included  in  the  evaluation 
of  the  thrust  and  the  -specific  impulse  of  the  engine . 
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Definition  of  Control  Volumes: 
ipdy/KL 


P  I  P  +  dP 

- v  '  j  v  +  dv 

P  I  (l)  j  a  +  dp 
P  |  din  P  +  dp 

* - - 

— *—  v  v  (2)  I  v  +  dv  ' 


I 

|_ _ J 


vN  Z?' 

\®17L- 


i is  ,  vt 

N  1 


2  *  ,,2 


Sectional  Vieiv 


Equations  for  control  volume (T) 
Continuity  Equation  in  (l) 


pvAc  +  dy  =  (p  +  dp)(v  +  dv)Ac  + 


d(pvAc)  =  dy  -  ^2- 


+  ^2. 
c  N 


Continuity  of  the  Unburned  Species 


m© 


/  ,  din 

(pvApJp  +  -*-  dy  =  (p  +  dp)(v  +  dv)A(,(p  +  dp)  +  — 


/  \  “p  din 

d(PVAcp)  =  dy  .  —  ^ 


Combination  with  the  continuity  equation  yields 


pvAcdp 


=  (1  -  P) % 


y- Component  of  Momentum  Equation  in 


I 

_pv2ac  ■  S  Vv  dy  +  v  ^  +  (p  +  dP)<v  +  dv)2  Ac 


PAe  -  (P  +  dP)  V 


d(pv2Ac)  +  Ac  dP  =  ^  Vw  dy  -  v 

c  ^  NL  N 


Combination  of  this  with  the  continuity  equation  gives, 
'pv  Ac  dv  +  Ac  dP  =  —2-  (Vv  -  v)  dy 


Energy  Equation  in  (l) 


pv  Ac(h  +  i.  v2)  +  (hp  +  i.  Vw2 )  dy 
2  NL  2 


(p  +  dp)(v  +  dv)Ac[h  +  dh  +  i-  (v  +  dv)2  ]  +  (h  +  iv2) 


d[pvAc(h  +  i  v2)  ]  =  -£  (hp  +  i  V?)  -  ^.(h  +  i-  v2) 


Combining  this  with  the  continuity  equation  results  in 


pvA,.  d(h  +  |  v2)  =  ^  (hp  +  |  Vw2  -  h  -  |  v2)  dy 


For  a  mixture  of  unburned  and  burned  propellant  the  enthalpy  is  given  by 


h  =  phA  +  (1  -  n)hg 

oT 


=  H 


CpA  dT  +  hA 


(o) 


Lo 


+  (1  -  n) 


J\  « *  V°’ 

-o 


C  T  +  MhA(o)  +  (1  -  u)hB(0) 


IT 


since  it  has  been  assumed  that  Cp^ 


Cp  ■  constant . 


=  CPB 

Now  Q  is  defined  by  the  relation, 


h 

hA 

Therefore ,  the  enthalpy  is  given  by 


A  *  $+*B(0) 


h  =  Cp  T  +  pQ  +  hB 


and  the  energy  equation  can  be  written  as 


(o) 


pvAcd(CpT  +  pQ  +  hB(o)  +  -/  )  =  (CpTp  +  Q  +  hB(0)  +  l  Vv2-CpT 


NL  '  P  P 
-  pQ  -  hB(0)  -  |  v2) 


Combination  with  the  continuity  of  species  relation  gives 


OvAcd(0pT  ♦  |  ,*)  ,  (0pTp  ♦  |  Vv2  -  0PT  -  i  ,.) 


Equations  for  control  volume  (?)  : 
Continuity  Equation  in  (2 ) 


J  J pvdAg  +  ^  =  J J  ( P  +  dp)(v  +  dv)dAg  + 
As  As 


*//"“*  -  £  -  £ 
As 


Since  it  is  assumed  that  the  fluid  flow  in  control  volume (5) is  quasi -one 
dimensional  in  the  axial  direction  it  is  clear  that  the  following  inequality 
must  be  satisfied, 


II 


PvdA 


[As 


dmt/N 


«  1 
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This  condition  can  he  checked  when  a  solution  to  the  equations  i6  obtained. 
Then  the  continuity  equation  becomes 

din  *  dmt 

Energy  Equation  ia  © 


ff  Pv(h+|v2)dA8+^(h+|v2) 


As 


P  +  dp)(v  +  dv)(h  +dh  +  |  (V  +  dV^ldAg  +  ^  (ht  +  J  VtC) 


A« 


d  J  J  Pv<h  +  I  Va)dAg  -  ~  (b  +  |  ^  (ht  +  |vta) 


Again,  it  is  assumed  that 

d£/pv(h  +  |v2)dAB 


N 


(bt  +|vt2) 


«  1 


Obtaining  as  a  result 


TT  (h  +  lv‘> 


which  becomes 


(o  )  1  9  (o  )  1  9 

CpTt  +  HQ  +  hB  +|  (at2  +  v  )  *  CpT  +  mQ  +  hB  +  -  v 


or 


3  2  2 

a*  =*  - —  a 

^  7+1 


The  incremental  mass  flow  out  of  the  throat  is  given  by 

<&t 
N 


—  ”  P-tatit(iy 
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and  since  it  is  consistent  with  the  assumptions  to  regard  the  fluid  flow  in  the 
control  volume  as  being  isentropic, 


%  -  ' 


7  +  1 


dy 


Therefore ,  the  equations  governing  the  flow  in  control  volume  (T)  may  be 
written  as 

ULk  1 


d(pvAc)  = 


NL 


pa 


2(7  -  1) 


dy 


PVAjjdH  =  (1  .  p)  ^  dy 


pvAcdv  +  AcdP  =  ^  (Vw  -  v)  dy 
pvAcd(CpT  +  |  v2)  =  (CpTp  +|VV2  -  CpT  -  |  v2)dy 


p 

Ko 

=  p  —  T 
m 

CpT 

a 

”7-1 

a2 

=  TL 

In  non-dimensional  form  the  equations  become 


d(p?)  =  (G-^-I)dH 

Mq  a 

pvdp  =  (l  -  p)GdT) 
pvdv  +  jtfT  dP  =  G(Vv  -  v)di, 


pvd(T  +  1Lj±  M02  v2)  =  G(Tp  +  -  f .  -  Z-I_i  Mq2v2  )drj 
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with  the  definitions, 


T  =  a2 
a2  =  P/p 


G 


p  V  A„N 
o  c 


B 


7  +  1 

—2 -  2(7  -  1) 

(r  +  1)  Ac 


If  P,  v,  "£,  and  p  are  chosen  to  be  the  dependent  variables  one  can  write 
the  equations  in  the  form, 


0  +  '  0 


+ 


0 


/Pv\dp 

Va2  /dq 


(1  -  n)G 


G(Tp  +  2L^i  Mr)2Vw2  -  a2  -  ^MQV) 


The  determinant  of  the  coefficients  is  given  by 


I 


D 


0 


v 


0 

P 


0 

2Pv 


Pv 

a2 

0 


Pv 

a2 

(y  -  Dm/pv* 

a2 


0 

2Pv 

a 


0 

0 


D 


2P3v2 

yM02as 


(1  -  M2) 


To  complete  the  statement  of  the  problem  for  the  complete  mixing  case  It 
is  necessary  to  specify  the  end  conditions.  At  tj  =  0  by  definition, 

fo  =  \>  =  &o  =  ^o  =  0 

and  at  tj  =  1  the  jump  conditions  across  the  detonation  wave  give  the  relations, 

Vi  =  b  +  'jb2  +  c 

Pi  =  1  -  7M02(v1  -  1) 

Si  =  |  Cl  +  yMQ2(i  -  Vi  >]vi 


where  b  and  c  reduce  to. 


b 


i  +  yMo 

(y  +  DMo2 


c 


2 

(y  +  Dm,,2 


i] 


6.  Equations  for  the  Case  of  No  Mixing  of  the  Burned  and  Unburned  Propellant 

Here  it  is  assumed  that  the  burned  and  unburned  gases  are  separated  by  an 
interface  which  is  stationary  with  respect  to  the  moving  waves  and  across  which 
no  mass  or  energy  transfers.  The  burned  and  unburned  gases  are  assumed  to  have 
distinct  values  of  specific  heat  ratios  and  molecular  weights,  namely,  73,  7^, 
Sg  and  m^.  In  other  respects  the  treatment  follows  that  of  the  complete  mixing 
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case.  Control  volume  fl)  is  split  by  the  interface  into  two  sections,  A  and  B 
corresponding  to  the  unturned  and  the  burned  propellants,  respectively.  In 
each  section  the  flow  is  assumed  to  be  quasi-one  dimensional  in. the  circumerentiftl 
direction.  The  pressure  gradient  across  the  interface  is  taken  to  be  zero.  The 
slope  of  the  interface  is  assumed  to  be  small.  Since  the  slope  of  the  interface 
is  given  by 

dx  _  u 
dy  ~  V 

i.e.,  the  interface  must  be  a  streamline.  This  assumption  is  consistent  with  the 
previous  condition 

V-«* 

The  treatment  of  control  volume  (sT)  is  identical  to  that  of  the  complete 
mixing  case,  and  the'  results  can  be  stated  as 


di^  =  dn^ 


IT  *  Vb  (^Tt)2(7b  '  X)  **  dy 


7b 


?b  +  1 


Definition  of  control  volumes: 

n^dy/NL 

0  _ /  /  /  //  //  / 

/////✓/////// 

v.  * 

P+dP 

x+dxj  VA  +  dvA- 
PA  +  dPA 


interface 
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Equations  for  control  volume  (a): 
Continuity  Equation  in  (a) 


PAvA*cx  +  ^  dy  =  (pA  +  dpA)(vA  +  dvA^c^ 

m 

d(pAVcx)  = 


Upon  integration  this  becomes,  since  x  =  0  when  y  =  0, 


PAVAJCX  =  Ey 


y -Component  of  the  Momentum  Equation  in  (a) 

t 

-  PAVA%X  -  rn  Vvdy  +  (pA  +  dPA)(vA  +  dvA)2lc' 


=  Plcx  -  (P  +  dP)lc(x  +  dx)  +  Picdx 


d(PAvA2*cx)  +  '*cxdP  =  NL  Vv  dy 


Combination  with  the  continuity  equation  gives 


-E-  y  dv.  +  /  xdP  =  -2  (V„  -  vA)  dy 
NL  A  c  NL  W  A 


Energy  Equation  in  © 


pMx^a  +  2-vaz)  +^(hP  +  |C 


=  (pA  +  dpA)(vA  +  dvA)/c(x  +  dx) [hA  +  dhA  +  | 


d  fPAvAfcx(hA  +  \  VA2 )  1  =  £  (hp  +  \  Vv2 ) 


+  dx) 


x  +  dx) 


dy 

A  +  dvA)2] 
dy 
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Integration  gives 


PAvA/cx(hA  +  |  v/' 


(h  +  i  V  2)  y 
NL  P  2  v  1  7 


vhich,  "by  the  continuity  equation,  simplies  to 


CPATA  +  Q  +  hB 


(o)  ,  l 


(o) 


.  1  „  2 


?  VA  =  CPATp  +  Q  +  hB  +  o  Vw 


Va  *  k  Vs  -  Vp  *  J  Vv2 


Equations  for  control  volume  (iT): 
Continuity  Equation  in  (b) 


dm 

•  pBvB^c^xn  "  x)  =  ~  +  (PB  +  dpB)(vB  +  dvB)ic(xn  -  x  -  dx) 

dm_ 

aiPeVB^^  -  =  -  — 

y -Component  of  the  Momentum  Equation  in 

2  dmn  2 

-  PBVB  J'cK  ‘  x)  +  IT  VB  +  (pB  +  dPBHvB  +  dvB)  -f  (x^  -  x  -  dx) 


=  Pic(xn  -  x)  -  (P  +  dP)!^^  -  x  -  dx)  -  Pfc  dx 


dm 


4[PBvB2*c(xn  -  X)1  +  ic(*n  "  x)dP  “  "  ~  VB 


Combination  with  the  continuity  equation  gives, 

PBvBdvB  +  dp  -  0 

Energy  Equation  in 

1  *1  2 

PBVc(xn  -  x)fhB  +  -  vB  ]  =  — [hB  +  g  vB  ] 

+  (PS  +  dPBHvB  +  dvB)fc(xn  -  x  -  dx)  [hg  +  dhB  +  |  (vB  +  dvB)2] 
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d[pBvBic(xn  -  x)(hB  +  |  vB2)] 


TT(hB+lvBS) 


Substituting  in  the  continuity  equation  yields 

d(hB  +  \  VB^  =  ° 
which  may  be  integrated,  to  give 


CPBTb  +  V  +  2  VB 


|  VB0  +  C-T*-  +  h' 


(o) 


PB  Bo  B 


CPBTB  +  2  VB2 


CPBTB0  +  \  VBo2 


Therefore,  the  pertinent  equations  for  the  case  of  no  mixing  are 


ih 

PaVcX  = 


y  dvA  +  ^  ^  (Vv  -  vA)dy 

JJL  A  JML. 


V*  *  |  ’a*  ■  cp.tp*Jv. 


-  -  iy* 

PA  P  2  w 

7fi  +  1 


d[pBvBic(xn  -  x)]  =  '  PB®B  (^TTl)  2(7B  ’  X)  **  47 


pBvBdvB  +  dP  =  0 


CPB  TB  +  I  tb“ 


cpbtBo  *  I X2 


Ro  Bo 

P  -  p4^ta  =  PB  =-  tB 


mA 


®B 


cpata 


7k 


-  1 


C  T _ 

58  5  "  7B  -  1 
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7jl  ' 

pa  ’ 


*B 


V 

pb 


In  non-dimensional  form  the  equations  become 

¥a5  =  GT> 


G^a  +  ^  dP  =  °(VW  -  VA)d T, 


VA,ar^»*  ■  SA  *  !s—i  *  * 


.  V  2 
‘o  vw 


d[pBVB(l  -  |)  ] 


BP 

^o®B 


dq 


^B7BdvB  + 


0 


=  *B  -  1 


2—2 


Tb  +  2 - Mo  VB  =  1  + 


rB  -  1 


p  "  PATA  ^  “  %*B 


C«.T 


J’B 


-  1 


a,2 


PaaA  7a  -  1  A  ’ 


%  *  aB* 


7a  £ 
7B  PA 


•»a 


The  momentum  equation  for  control  volume  (IT)  can.  be  integrated  by  finding 
a  relation  for  Pb  in  terms  of  and  P. 


P 


27 


1 


Mo2(vb2  -  1) 


aB2 


7b  -  1 


Hence,  the  momentum  equation  becomes 


1 

2vBdvB 

+  1 

~&P  | 

(rB  -  Dm0» 

Vj2  -  1  -  2 

(7B  -  l)Mo2J 

ITJ 

Giving  upon  integration 


/ 


7b  -  1 


f  *  <*b  -  'K  t 


?  ^ 
2  ix  -  r 


which  leads  to  the  relations 


aB 


7B  -  1 

=  p  ^7B 


-1/7B 
P 

We  have  two  remaining  differential  equations: 


=  P 


G^A  +  ^dP  =  G(Vw-vA)dT, 


d[pBvB(l  -  |)] 


fj-dt, 

Ho  aB 


which  can  be  written  to  include  just  two  dependent  variables . 

Boundary  conditions: 
at  Tj  =  0, 
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where  in  this  case 
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It  is  implied  in  the  no-mixing  flow  model  that  at  r)  =  l,  5=1.  For  this 
condition  to  hold  the  continuity  equation  for  control  volume  (a)  and  the  first 
Jump  condition  give  G  =  1. 


7'.  Conclusions 

The  assumptions,  conditions,  and  resulting  equations  for  an  analytical 
model  of  the  rotating  detonation  wave  engine  have  been  presented.  The  validity 
of  this  model  rests  upon  the  premise  that  the  fluid  flow  is  quasi-one  dimensional 
when  vieved  in  the  system  of  the  rotating  detonation  waves .  The  effects  of  tur¬ 
bulent  mixing  of  the  burned  and  unburned  propellants  are  to  be  estimated  by  ap¬ 
plying  the  alternative  limiting  conditions  of  complete  mixing  and  no  mixing. 

In  both  cases  the  resulting  equations  are  a  system  of  a  first  order  ordinary 
differential  equations  with  initial  values  related  to  the  hydrodynamic  Jump  con¬ 
ditions  across  the  detonation  waves . 

Future  work  will  consist  of  obtaining  solutions  to  the  equations,  either 
exact  or  approximate,  correlating  the  theory  with  data  from  the  100-lb  engine 
tests,  and  using  these  results  to  aid  the  1000-lb  engine  design  effort. 
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II.  EXPERIMENTAL  STUDIES 


A.  THE  GASEOUS,  100- LB  THRUST  MOTOR 

The  revision  of  the  existing  100-lb  thrust  motor  began  with  modifications 
which  make  it  possible  to  change  the  injector  faces,  throat  sections,  and  the 
axial  length  of  the  combustion  chamber  quite  rapidly.  These  modifications 
will  make  it  possible  to  study  the  effect  of  L*,  the  characteristic  length, 
on  the  performance  of  the  100-lb  thrust  motor.  Also,  the  interchangeability 
of  the  injector  face  on  the  modified  100-lb  thrust  motor  will  enable  the  ef¬ 
fect  of  mixture  ratio  and  injector  pressure  to  be  determined. 

Originally  the  motor  was  located  in  a  pit  11  feet  deep  in  the  center  of 
the  isolation  cell.  However,  since  the  initiation  of  the  contract,  the  motor 
has  been  moved  to  ground  level.  The  relocation  of  the  motor  will  facilitate 
day  by  day  modifications  and  checks  that  must  be  carried  out,  also,  the  new 
location  should  make  it  possible  to  obtain  photographs  of  the  exhaust  jet  more 
easily. 

The  propellant  feed  system  has  also  been  modified  by  the  installation  of 
a  new  Hoke  pressure  regulator  which  appears  to  have  a  zero  drift  characteris¬ 
tic.  The  original  pressure  regulator  was  being  operated  above  its  maximum 
pressure  limit  which  resulted  in  an  extensive  drift  in  the  injector  pressure 
between  runs. 

The  fire  controller  has  also  been  extensively  modified.  With  the  instru¬ 
mentation  in  its  original  configuration  there  appeared  to  be  an  interconnec¬ 
tion  between  the  fire  controller  or  time-delay  spark  unit  and  the  oscilloscope 
trigger  circuit.  The  revision  of  the  fire  controller  and  the  relocation  of 
both  the  oscilloscope  and  the  time  delay  spark  unit  have  apparently  eliminated 
this  interconnection. 

The  original  instrumentation  has  been  reassembled  and  relocated.  This 
entire  system,  solenoid  valves,  injector  pressure  transducers  and  associated 
instrumentation,  oscillograph,  fire  controller,  time-delay  spark  unit,  etc., 
has  been  checked  by  sequencing  the  motor  with  gaseous  nitrogen,  see  Fig.  4, 
a  schematic  diagram  of  the  revised  test  setup.  Also  see  Figs.  5  and  6  for 
photographic  details  of  the  test  setup.  One  of  the  oscillograph  records  for 
this  shakedown  phase  is  shown  in  Fig.  7<  The  record  shows  the  essential  mech¬ 
anical  operations  plus  the  injector  pressure  and  the  ignition  of  the  motor. 

The  Kistler  pressure  transducer  instrumentation  used  to  record  the  chamber 
pressure  as  a  function  of  time,  was  installed  in  the  100- lb  thrust  motor  and  a 
live  run  made.  Since  these  pressure  traces  were  completely  unsatisfactory,  the 
chamber  pressure  instrumentation  has  been  removed  from  the  motor  and  installed 
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in  a  detonation  tube.  It  is  believed  that  this  will  facilitate  the  detection 
of  the  source  of  a  high  frequency  oscillation,  that  has  been  superimposed  on 
the  pressure  trace,  more  efficiently  than  could  be  done  in  the  100-lb  thrust 
motor,  This  configuration  will  also  make  it  possible  to  obtain  dynamic  pres¬ 
sure  calibrations  of  the  transducer  by  comparing  the  experimental  results  in 
the  tube  with  the  theoretical  results  predicted  from  non- steady,  one-dimen- 
sional  flow  theory.  An  example  of  such  a  pressure  trace  utilizing  the  Kistler 
pressure  transducer  with  the  voltage  displayed  on  the  Tektronix  Model  555  os¬ 
cilloscope  is  shown  in  Fig.  8.  The  horizontal  sweep  time  is  50  u-sec/cm  (cm  = 
large  division) .  While  it  is  apparent  that  some  high-frequency  oscillations 
in  the  trace  occur  that  are  most  pronounced  immediately  after  the  passage  of 
the  first  wave,  the  vertical  deflection  (pressure)  decreases  to  about-  one- 
third  its  value  immediately  behind  the  wave  in  about  400  p-sec.  This. is  in 
good  agreement  with  theory.  The  second  pressure  pulse  (at  about  460  p-sec) 
is  a  reflected  shock  wave  from  the  diaphragm  located  at  the  end  of  the  tube. 

The  thermocouples  used  to  measure  the  wall  temperature  of  the  combustion 
chamber  will  probably  be  checked  and  dynamically  calibrated  in  much  the  same 
way  as  the  Kistler  pressure  transducer. 

During  the  period  in  which  the  pressure  transducer  is  being  evaluated 
and  calibrated,  two  high-speed  framing  cameras  and  one  time  exposure  camera 
are  being  installed  in  the  isolation  cell.  These  high-speed  framing  cameras, 
a  16  mm  Fastax  and  a  l6  mm  Beckman  and  Whitley,  will  obtain  photographs  of 
the  exhaust.  The  velocity  of  the  detonation  wave  and  also  any  acceleration, 
whether  a  positive  or  negative,  can  be  determined  from  these  photographs. 

The  photographs  from  the  time  exposure  camera  will  ascertain  the  shape  of 
the  Jet  and  help  shed  some  light  on  the  velocity  pattern  of  the  exhaust  Jet. 

These  cameras  have  been  installed  in  the  test  cell,  and  the  electronic 
unit  necessary  to  sequence  them  is  nearing  completion. 


Special  Instrumentation 

(a)  Chamber  Pressure .—The  chamber  pressure  instrumentation  consists 
of  a  Kistler  Model  605  quartz  crystal  pressure  transducer,  the  output  of 
Which  is  fed  into  a  Kistler  Model  566M  charge  amplifier.  The  output  of  the 
amplifier  is  applied  to  the  CRT  of  a  Tektronix  Model  555/  dual-beam  oscillo¬ 
scope  and  photographed  with  a  Polaroid  oscilloscope  camera.  The  entire  sys¬ 
tem  should  have  a  rise  time  <  2  p-sec. 

(b)  Combustion  Chamber  Wall  Temperature .  —The  wall  temperature  of  the 
combustion  chamber  will  probably  be  measured  with  a  Nanmac  film  thermocouple 
(10  p-sec  rise  time)  in  conjunction  with  a  voltage  amplifier  and  the  Tek¬ 
tronix  Model  555  oscilloscope.  The  fast  rise  of  this  type  of  thermocouple 
will  make  it  possible  to  predict  steady  state  heat  transfer  rates  with  a  high 
degree  of  accuracy. 
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(c)  Exhaust  Jet . —The  Fastax  l6  mm  framing  camera  is  rated  at  approximately 
7000  frames  per  second  which  will  give  about  one  frame  per  revolution  of  the 
detonation  wave.  This  camera  also  has  sufficient  film  capacity  (100  ft)  to  re¬ 
cord  the  entire  run  for  run  times  of  less  than  one  second. 

The  Beckman  and  Whitley  high-speed  framing  camera  has  a  framing  rate  of 
about  26,000  frames  per  second  but  due  to  its  limited  film  capacity  (33-7/8 
in.),  only  approximately  8  milliseconds  of  the  run  can  be  observed.  Even 
though  only  a  short  portion  of  the  run  can  be  recorded, the  high  framing  rate 
will  yield  much  more  accurate  data  on  wave  velocities  and  accelerations  than 
is  possible  with  the  Fastax  camera. 

The  photographs  from  the  time-exposure  camera  will  show  the  shape  of  the 
jet  which  in  turn  should  provide  some  information  concerning  the  nature  of  the 
expansion  processes  of  the  exhaust  gases. 


B.  TEMPERATURE  AND  PRESSURE  EFFECTS  ON  HYDROGEN- OXYGEN  DETONATION  VELOCITIES 

The  experimental  setup  has  been  completed.  The  additional  items  of 
equipment  inlcude  the  following: 

1.  A  stainless  steel,  vacuum- jacketed  vessel  to  contain  the  coiled 
detonation  tube  and  the  low  temperature  liquid  (see  Fig.  9). 

2.  An  exhaust  fan  to  eliminate  noxious  fumes  and  combustible  gases 
from  the  isolation  cell, 

3.  A  CMC  Model  727-BN  time  interval  counter  capable  of  recording 
time  intervals  to  an  accuracy  of  +  0.1  n-sec. 

It  was  decided  that  three  different  cooling  fluids  would  be  utilized. 

The  first,  normal  propyl  alcohol,  used  by  Moyle1 5  in  his  experiments,  would 
be  utilized  from  room  temperature  down  to  near  its  freezing  point,  about 
l60°K,  The  second,  isopentane  wo'xld  be  utilized  tc  cover  the  range  in  tem¬ 
perature  from  160 °K  to  about  10r.°K,  where  it  freezes.  The  range  from  102“K 
down  to  90*K  (the  vapor  saturation  temperature  of  oxygen  at  one  atmosphere) 
would  be  covered  utilizing  liquid  nitrogen  pressure-regulated  at  a  particular 
pressure  to  obtain  the  temper it-ure  corresponding  to  that  vapor  pressure  of 
nitrogen.  It  is  to  be  noted  that  cold  nitrogen  vapor  initially  at  about  80“K 
is  bubbled  through  the  bath  to  cool  either  of  the  two  fluids,  normal  propyl 
alcohol  or  isopentane  in  the  higher  temperature  range.  Table  VI  shows  the 
temperature  and  pressure  regions  to  be  covered  in  the  experiments  and  the  cool¬ 
ing  fluid  to  be  used  in  the  various  temperature  regions. 

The  stainless-steel  coiled  detonation  tube  utilized  in  the  experiments 
has  an  I.D.  of  0.25  in.  similar  to  the  one  used  by  Moyle.  It  is  approximately 
25  ft  long  with  the  middle  12  ft  submerged  in  the  bath  and  used  for  the  vel- 
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ocity  measurement.  Three  ionization  probes  are  installed  6  feet  apart  allowing 
the  detonation  velocity  to  be  measured  over  two  6-foot  intervals  in  the  same 
run,  affording  a  check  that  the  wave  has  reached  its  steady- state  Chapman- Jou- 
guet  velocity. 

Preliminary  tests  were  made  in  the  coiled  tube  to  check  the  fact  already 
noted  by  Moyle,  i.e.,  that  the  detonation  velocity  measured  in  a  0.25  in.  I.D. 
tube  coiled  to  a  ten- inch  diameter  is  very  close  (about  0.5$  less)  than  the 
detonation  velocity  measured  in  a  straight  tube.  These  were  also  run  to  es¬ 
tablish  that  at  room  temperature,  H2-02  detonations  propagated  at  a  constant 
Chapman-Jouguet  velocity  over. the  12-foot  length  of  the  coil  situated  inside 
the  bath. 

The  first  series  of  experiments  were  performed  to  check  the  results  of 
Moyle  at  one  atmosphere  and  at  different  mixture  ratios  of  hydrogen  and  oxy¬ 
gen  from  room  temperature  down  to  l60°K.  The  experiments  were  performed  in 
the  same  manner  as  Moyle:  normal  propyl  alcohol  was  utilized  as  the  cooling 
fluid  with  cold  nitrogen  vapor  used  to  cool  the  alcohol  to  the  desired  tempera¬ 
ture.  After  each  run  the  coil  was  removed  from  the  bath  and  allowed  to  warm 
up  to  room  temperature  to  melt  the  combustion  products  (HaO)  that  had  frozen 
in  the  tube.  The  tube  was  dried  under  a  vacuum  and  recharged  with  the  H2-02 
mixture  and  then  replaced  in  the  bath.  The  results  of  these  experiments  are 
shown  in  Fig.  10.  Although  the  same  mixture  ratios  were  not  employed,  it  is 
believed  that  the  experimental  verifications  of  Moyle's  data  can  now  be  es¬ 
tablished.  .  It  had  been  hoped  that  in  order  to  make  a  large  number  of  runs 
in  a  minimum  of  time,  the  coiled  tube  could  remain  in  the  cooling  bath  be¬ 
tween  runs  and  the  frozen  combustion  products  eliminated  by  one  of  the  follow¬ 
ing  methods: 

1.  Purge  coiled  tube  with  dry  N2  immediately  after  run. 

2.  Sublime  the  ice  by  means  of  a  vacuum. 

3.  Utilize  a  helium  shock- tube  driver  section  upstream  of  the  coil 
and  initiate  the  detonation  wave  with  shoc.t  wave  with  the  re¬ 
sulting  helium  piston  gases  purging  the  tube  of  HgO. 

The  first  two  methods  appeared  to  have  no  success.  Elimination  of  the 
ice  was  deemed  successful  if  a  series  of  runs  made  at  one  low  temperature 
yielded  the  same  detonation  velocity.  Since  the  velocity  was  observed  to 
change  by  about  5$  it  was  concluded  that  ice  was  still  accumulating  on  the 
inside  walls  of  the  tube.  The  third  method,  employing  the  helium  shock  tube 
driver  appeared  to  be  successful  inasmuch  as  the  detonation  velocities  were 
the  same  for  several  runs  at  moderate  sub-freezing  temperatures  and  the  re¬ 
sults  agreed  with  the  control  experiments  performed  by  removing  the  tube  from 
the  bath  for  each  run.  However,  at  lower  temperatures  it  appeared  as  though 
the  HgO  was  not  being  removed  rapidly  enough  by  the  He  piston  gases  and  ice 
was  beginning  to  form. 

The  following  decision  has  been  made:  the  tedious  method  of  removing  the 
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coil  from  the  bath  will  be  utilized  and  the  pressure  and  temperature  region  shown 
in  Table  VI  will  be  investigated  down  to  the  freezing  point  of  isopentane  (112°K). 
Since  the  bath  will  be  at  one  atmosphere  for  these  runs,  the  cover  will  not  have 
to  be  unbolted  and  rebolted  after  every  run,  a  time  consuming  feature  if  LN2  is 
employed  as  the  coolant.  By  utilizing  this  technique,  it  will  not  be  possible 
to  check  the  region  from  90°K  to  112°K.  However  only  a  very  small  amount  of 
the  data  was  to  be  recorded  in  this  range  anyway  (as  is  evident  by  examining 
Table  Vi),  due  to  the  higher  saturation  temperatures  of  oxygen  at  increased 
pressure. 


C.  DETONATION  THROUGH  HETEROGENEOUS,  LIQUID-GAS  MEDIA 

A  system  involving  a  condensation  chamber  has  been  fabricated  and  instru¬ 
mented  and  will  be  used  to  appraise  the  characteristics  and  limitations  of  the 
condensation  system  for  droplet  formation.  This  system  will  utilize  condensa¬ 
tion  of  one  component  of  a  mixture  due  to  the  temperature  decrease  associated 
with  an  isentropic  expansion  of  the  gases  in  the  chamber.  At  the  same  time  it 
will  provide  a  means  for  determining  the  capabilities  of  the  schlieren  system, 
which  will  be  used  on  the  detonation  bube  as  well.  Figure  11  shows  a  view  of 
the  condensation  chamber  with  glass  windows  as  well  as  one  leg  of  the  schlieren 
system  involving  6  in.  diameter  optics. 

The  schlieren  system  is  almost  completed  and  is  being  tested  at  various 
degrees  of  magnification  using  supplemental  lenses.  It  is  believed  that  the 
well  known  limitation  of  film  speed  and  resolution  will  limit  the  magnifica¬ 
tion  to  between  five  and  ten  times  natural  size  utilizing  for  exposure  the 
light  that  is  available  from  the  high-voltage  (25  kilivolt)  short  duration 
(about  0.1  n-sec)  light  source. 

The  detonation  tube  has  been  fabricated  in  its  component  pieces  and  is 
ready  for  final  assembly.  Tests  on  mixtures  containing  the  very  small  drop¬ 
lets  produced  by  the  condensation  technique  will  be  performed  first.  The  in¬ 
strumentation  involved  for  this  phase  will  consist  of  the  schlieren  system 
utilizing  spark  photography,  ionization  probes  to  determine  detonation  wave 
velocity,  and  a  thermocouple  and  pressure  transducer  to  determine  the  state 
of  the  mixture  prior  to  detonation. 

An  injection  system  capable  of  producing  larger  droplets  of  either  fuel 
or  oxidizer  is  under  consideration. 


D.  GEOMETRICAL  TESTS 

The  experimental  setup  described  in  the  last  report^  has  been  completed. 
The  following  changes  have  been  made: 
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1.  Glow- plug  ignition  ie  utilized  instead  of  spark  ignition  to 
eliminate  the  spark-discharge  interference  with  the  electronic 
equipment . 

2.  Plexiglas  windows  with  inferior  optical  quality  were  utilized 
initially  in  the  test  section  due  to  the  difficulty  Of  fabri¬ 
cating  curved  glass  windows. 

3.  A  conventional  schlieren  system,  employing  8-inch  diameter 
mirrors  with  a  short  duration  (0.1  u-sec)  high  voltage  (25 
kilovolt)  capacitor-discharge  light  source  was  employed . 

4.  A  thyratron-triggered  time  delay  unit  consisting  of  a  variable 
R-C  time  delay  circuit  is  utilized  to  trigger  the  light  source 
at  the  proper  time. 

5.  A  CMC  Model  757-BN  time-interval  counter  is  used  to  measure 
the  wave  velocity  in  the  straight  or  curved  tube. 

Figure  12  shows  a  view  of  experimental  apparatus  and  Fig.  13  shows  a  detailed 
photograph  of  the  curved-tube  section. 

Preliminary  experiments  have  been  performed  utilizing  the  curved  test 
section.  The  internal  dimensions  of  the  rectangular  channel  are  0.5  in.  in 
the  radial  direction  and  0.375  in.  vide.  The  curved  section  has  a  7-5  in. 
diameter,  similar  to  the  100-lb  thrust  motor. 

The  experimental  procedure  is  as  follows t 

1,  The  desired  mixture  of  Hg  and  02  is  prepared  in  the  mixing  sphere 
and  allowed  to  stand  for  24  hours  to  insure  uniform  composition. 

2.  The  detonation  velocity  is  measured  in  the  straight  tube  to  as¬ 
certain  the  exact  mixture  ratio  (utilizing  the  extensive  data 
by  Moyle ) . 

3,  The  straight  tube  with  the  curved  section  is  charged  with  the 
desired  mixture  and  ignited. 

4.  The  wave  triggers  the  time-delay  unit  which  in  turn  triggers 
the  light  source  taking  a  schlieren  photograph  of  the  wave  in 
the  curved  section  of  the  tube. 

Figures  l4a  and  l4b  are  Bpark-schlieren  photographs  taken  of  a  Hg-0a  detona¬ 
tion  with  a  mixture  ratio  of  60^  H2  by  volume.  The  photographs  were  obtained 
on  Polaroid  Land  film  and  two  different  magnifications  were  employed.  The 
wave  in  each  case  has  passed  through  about  3/4  of  the  curved  section.  Several 
interesting  phenomena  can  be  observed  on  the  photographs.  The  leading  edge  of 
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the  wave  front  appears  to  he  at  a  slight  angle  to  the  radial  direction,  possibly 
showing  the  effect  of  the  greater  distance  the  wave  must  travel  at  the  outer 
radius.  Also,  a  multiply-reflected  wave  is  visible  in  the  burned  gases  behind 
the  initial  wave  front.  This  effect  has  been  noted  in  earlier  experiments  per¬ 
formed  at  this  laboratory,  but  in  these  cases  the  reflected  wave  in  straight 
tubes  were  always  noted  during  the  starting  transient  of  the  detonation,  i.e., 
shortly  after  the  transition  from  deflagration  to  detonation.  These  disturb¬ 
ances  always  died  out  and  were  not  apparent  at  all  if  a  sufficient  distance 
was  allowed  for  the  wave  to  travel  after  its  transition  to  a  detonation  wave. 
Since  the  distance  between  the  point  where  the  wave  is  initiated  in  the  straight 
section  and  where  the  photographs  are  taken  in  the  curved  section  is  more  than 
adequate  for  this  disturbance  due  to  the  initiation  of  the  wave  to  dissipate, 
it  is  concluded  that  the  observed  disturbance  in  the  curved  section  is  caused 
by  the  effect  of  the  curved  channel  on  the  detonation  process.  Whether  this 
disturbance  dissipates  with  distance  or  is  continually  being  generated  must 
be  established.  It  is  possible  that  this  disturbance  could  have  an  effect  on 
the  detonation  wave  structure  and  propagation  velocity,  if  it  persists.  Uti¬ 
lizing  many  photographs  of  different  waves  at  varying  distances  along  the  curved 
tube  will  check  this  effect  in  addition  to  establishing  the  effect  of  curvature 
on  wave  velocity  and  detonation  limits. 


III.  STUDY  PLANS  FOR  THE  NEXT  QUARTER 


The  theoretical  studies  described  in  Sections  I-A  and  I-B  of  this  report 
will  be  terminated  early  in  this  period.  However  experimental  results  will  be 
utilized  later  to  check  both  theoretical  models.  The  study  of  the  theoretical 
model  describing  the  gas-dynamics  involved  with  the  rotating  detonation  wave 
engine  will  be  continued.  Approximate  solutions  of  the  equations  for  the  com¬ 
plete  mixing  case  as  well  as  for  the  case  without  mixing  will  be  attempted  at 
an  early  date.  This  will  be  done  to  aid  in  understanding  .the  experimental  data 
from  the  100- lb  thrust  motor  in  addition  to  providing  a  starting  point  for  more 
accurate,  detailed  numerical  solutions  that  could  be  obtained  at  a  later  date 
utilizing  a  digital  computer. 

A  complete  parametric  study  of  the  effects  of  injector  geometry,  pro¬ 
pellant  mass  flow  and  mixture  ratio,  and  chamber  axial  length  and  nozzle  throat 
size,  on  the  measurable  detonation  properties,  i.e.,  wave  velocity, chamber 
pressure,  and  wall  temperature  will  be  made  for  the  100-lb  thrust  motor.  These 
results  will  be  checked  with  the  theoretical  results  predicted  by  the  study 
described  in  Section  I-G  of  this  report. and  in  Reference  1. 

The  maximum  possible  range  of  pressure  and  temperature  effects  on  H2-O2 
detonation  velocities  will  be  studied. 

Further  experiments  will  be  carried  out  on  the  effects  of  curvature  and 
relief  on  detonation  wave  characteristics  in  the  curved  tube. 

Experiments  will  be  made  to  obtain  phenomenalogical  data  on  detonation 
waves  passing  through  heterogeneous  liquid-gas  media. 

Preliminary  and  then  final  design  will  begin  on  the  1000-lb  thrust  ro¬ 
tating  detonation  wave  engine  to  be  tested  at  the  Willow  Run  rocket  test  site. 
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TABLE  I 


THEORETICAL  AND  EXPERIMENTAL  CORRELATION  OF  CRITICAL  SIZE  OF  DROPLETS 
OF  DIFFERENT  FLUIDS  FOR  BAG-TYPE  DROPLET  BREAKUP 


Liquid 

Surface 

Tension, 

dyne/cm 

S 

Experimental 

Critical 

Velocity, 

ft/sec 

U 

Experimental 
Criti cal 
Size, 
micron 
(r)cr 

Calculated  Critical  Size 
by  Using  Eq.  (3)  and. 
Values  of  Experimental 
Critical  Velocity, 
micron 
(r)cr 

Burning  RP-1 

11 

92 

100 

75 

< 

Burning  RP-1 

11 

85 

100 

87.5 

§ 

Burning  RP-1 

11 

82 

100 

98.5 

o 

6 

Non-Burning  RP-1 

38 

95 

100 

240 

_Non-Burning  RP  1 

38 

125 

100 

l4o 

Dow  Corning  200-A  Fluid 

21.14 

73-5 

213 

246 

Dow  Corning  200-A  Fluid 

21.14 

93-0 

143 

156 

Dow  Corning  200-A  Fluid 

21.14 

121.3 

90 

103 

Dow  Coming  200-A  Fluid 

21.14 

148.3 

59 

60.5 

g 

Water 

71-97 

121.3 

200 

321 

5 

Water 

71.97 

198.8 

90 

119 

Water 

71.97 

282.0 

53 

59.5 

Methyl  Alcohol 

22.2 

74.4 

235 

263 

Jfethyl  Alcohol 

22.2 

121 

93 

101 

Water 

72.75 

84.3 

300 

665 

Water 

72.75 

109.5 

205 

377 

Water 

72.75 

157-3 

135 

180 

Water 

72.75 

238.5 

60 

73.5 

§• 

o 

u 

Methyl  Alcohol 

22.6 

60 

317 

415 

o 

Methyl  Alcohol 

22.6 

84.3 

165 

208 

Methyl  Alcohol 

22.6 

109.5 

125 

120 

Methol  Alcohol 

22.6 

157.3 

99 

55 

Experimental  data  for  group  (A) 
Experimental  data  for  group  (B) 
Experimental  data  for  group  (C) 


are  taken  from  Ref. 
are  taken  from  Ref. 
are  taken  from  Ref. 
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11. 

4. 
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TABLE  II 


THEORETICAL  AND  EXPERIMENTAL  CORRELATION  OF  CRITICAL  SIZE  OF 
BURNING  AND  NON-BURNING  DROPLETS  OF  RP-1  FOR  SHEAR-TYPE 
DROPLET  BREAKUP 


Liquid 

Surface 

Tension, 

dyne/cm 

S 

Experimental 
Critical 
Velocity, 
ft/ sec 

U 

Experimental 

Critical 

Size, 

micron 

(r)cr 

Calculated 

Critical  Size 
by  Using  Eq.  (5), 
micron 
(r)cr 

Burning  RP-1 

11 

83 

500 

152 

Burning  RP-1 

11 

55 

500 

3^5 

Non-Burning  RP-1 

38 

93 

500 

i+ll 

Non- Burning  RP-1 

38 

8o 

500 

550 

Experimental  data 

are  taken 

from  Ref.  6. 

TABLE  III 

THEORETICAL  SHEAR-TYPE  DROPLET  BREAKUP  TIMES  FOR  BURNING 
AND  NON-BURNING  DROPLETS  OF  RP-1  USING  EXPERIMENTAL  VALUES  OF 
CRITICAL  VELOCITY  AND  CRITICAL  SIZE 


Liquid 

Surface 
Tension, 
dyne/ cm 

S 

Experimental 
Critical 
Velocity, 
ft/ sec 

U 

Experimental 

Critical 

Size, 

micron 

(r)cr 

Breakup  Times 
Calculated  by 
Using  Eq.  (4-a), 
sec 

Burning  RP-1 

11 

83 

500 

0.1+89  x  10" 3 

Burning  RP-1 

11 

55 

500 

1.16  x  10" 5 

Non-Burning  RP-1 

38 

93 

500 

0.925  x  10-3 

Non-Burning  RP-1 

38 

80 

500 

1.29  x  10-3 

Experimental  data  are  taken  from  Ref.  6. 


1+2 


TABLE  IV 


THEORETICAL  SHEAR-TYPE  DROPLET  BREAKUP  TIMES 
OF  LIQUID  OXYGEN  DROPLETS  IN  A  STOICHIOMETRIC  H2-02  DETONATION 


Radius  of  the 
Droplets, 
r(microns) 

Reynold' s 
Number, 

Re 

Breakup  Time  Calculated 
by  Using  Eq.  (4-a), 
tb(sec) 

500 

30.8 

0.89  x  10" 7 

600 

37 

0.99  x  10"  7 

700 

43.1 

1.0  x  10"  7 

800 

49.2 

1.15  x  10"  1 

900 

55-6 

1.22  x  10" 7 

1000 

6l.6 

1.25  x  10“ 7 

Surface  tension  of  liquid  oxygen 
Detonation  velocity 
Mole-fraction  of  hydrogen 
Viscosity  of  gas  mixture 
Density  of  gas  mixture 
Density  of  liquid  oxygen 


13.2  dyne/cm  (at  T  =  90°K) 
9400  ft/sec 

O.667 

0.894  x  10-5  ”L 

5.2  x  10-i+  gm/cm^J 

0.4299  gm/ cm3  (at  T  =  154. 3°K) 


♦Corresponding  to  the  condition  behind  the  detached  shock. 


TABLE  V 

THE  VARIATION  OF  DETONATION  PARAMETERS  AND  HEAT 
FLUX  WITH  MIXTURE  RATIO 


X«2 

vD, 

ft/sec 

°R 

‘Iw 

Dissociation 
Neglected, 
Btu/ft^  sec 

8w 

Dissociation 
Included, 
Btu/ft2  sec 

0.40 

6,879 

5893 

1495  (^)1/5 

1902  (i)l/5 

0.50 

7,672 

6314 

1819 

2589 

0.60 

8,620 

6510 

2135 

3298 

0.76 

10,758 

6484 

2723 

3776 

0.80 

11,262 

6267 

2683 

3337 

=  mole-fraction  of  H2 
Vjj  =  detonation  velocity 
Te  =  temperature  behind  the  detonation 
qw  =  heat  flux  at  the  wall 
x  =  distance  behind  detonation  (ft) 
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Fig.  2.  Heat  flux  to  the  wall  of  the  rotating  detonation  wave  engine  versus  mixture 
ratio  with  and  without  dissociative  recombination. 


TEMPERATURE 


Fig.  3-  H2-O2  Chapman- Jouguet  detonation  and  combustion  ehambei 

temperatures  as  a  function  of  mixture  ratio. 
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Fig.  5-  Photograph  of  valve  and  injector  manifold  arrangement  for 
100-lb  thrust  motor. 


Fig.  6.  Photograph  of  downstream  end  of  100-lb  thrust  motor. 
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Fig.  7.  Oscillograph  record  of  essential  mechanical  operations,  injector  pressure 
and  ignition  of  100-lb  thrust  motor. 


Fig.  8.  Pressure  trace  of  H2-O2  detonation  wave  in  straight  tube  using 
the  Kistler  Model  603  pressure  transducer. 
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Fig.  12.  Photograph  of  experimental  equipment  associated  with  the  geometrical  effect  studies. 
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Fig.  13.  Photograph  of  curved  detonation  tube  section. 


Jig.  14.  Spark-schlieren  photographs  at  different  magnifications  of  two  6of,  Ha -02  detona¬ 
tion  waves  passing  through  the  curved  test  section. 
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